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Abstract 

Background: Biofeedback of body motion can serve as a balance aid and rehabilitation tool. To date, mathematical 
models considering the integration of biofeedback into postural control have represented this integration as a 
sensory addition and limited their application to a single degree-of-freedom representation of the body. This study 
has two objectives: 1) to develop a scalable method for incorporating biofeedback into postural control that is 
independent of the model's degrees of freedom, how it handles sensory integration, and the modeling of its 
postural controller; and 2) to validate this new model using multidirectional perturbation experimental results. 

Methods: Biofeedback was modeled as an additional torque to the postural controller torque. For validation, this 
biofeedback modeling approach was applied to a vibrotactile biofeedback device and incorporated into a two-link 
multibody model with full-state-feedback control that represents the dynamics of bipedal stance. Average response 
trajectories of body sway and center of pressure (COP) to multidirectional surface perturbations of subjects with 
vestibular deficits were used for model parameterization and validation in multiple perturbation directions and for 
multiple display resolutions. The quality of fit was quantified using average error and cross-correlation values. 

Results: The mean of the average errors across all tactor configurations and perturbations was 0.24° for body sway 
and 0.39 cm for COP. The mean of the cross-correlation value was 0.97 for both body sway and COP. 

Conclusions: The biofeedback model developed in this study is capable of capturing experimental response 
trajectory shapes with low average errors and high cross-correlation values in both the anterior-posterior and 
medial-lateral directions for all perturbation directions and spatial resolution display configurations considered. The 
results validate that biofeedback can be modeled as an additional torque to the postural controller without a need 
for sensory reweighting. This novel approach is scalable and applicable to a wide range of movement conditions 
within the fields of balance and balance rehabilitation. The model confirms experimental results that increased 
display resolution does not necessarily lead to reduced body sway. To our knowledge, this is the first theoretical 
confirmation that a spatial display resolution of 180° can be as effective as a spatial resolution of 22.5°. 

Keywords: Biofeedback, Sensory augmentation, Sensory substitution, Mathematical modeling, Postural control, 
Vibrotactile feedback, Multidirectional perturbations, Vestibular, Balance 



Background 

Biofeedback can be used to supplement or replace missing 
sensory information by providing the user with informa- 
tion via a functioning sensory modality. To date, indivi- 
duals with vestibular deficits and older adults have 
demonstrated improved balance when using electrotactile 
[1-6], vibrotactile [7-19], auditory [20-25], or multimodal 
[26,27] feedback displays of body motion during quiet or 
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perturbed standing and gait tasks. For example, vibrotactile 
feedback of torso tilt has experimentally shown reductions 
in the root-mean-square sway in subjects with unilateral 
and bilateral vestibular loss during both linear and rota- 
tional single-axis [10,11] and multidirectional [14,28] per- 
turbations of stance. Research is now underway to evaluate 
biofeedback, also referred to as sensory substitution or sen- 
sory augmentation, as a real-time balance aid and as a tool 
for balance rehabilitation. 

Biofeedback-related research aimed at improving balance 
has primarily used experimental rather than mathematical 
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methods to study the relationship between non-native 
feedback channels and postural control As demonstrated 
by numerous physiological studies [29-38], mathematical 
models can complement experimental work by allowing 
for design evaluation and optimization prior to human 
subject testing, explaining experimental findings and iden- 
tifying dominant underlying physiological mechanisms. 

The literature discusses different mathematical models 
of human posture with varying levels of complexity. The 
simplest model used to describe bipedal postural control 
is an inverted pendulum with feedback control [29-31], 
which is a one-link representation capturing a single de- 
gree of freedom. For example, Peterka used this model 
to explain the experimental differences in the sensori- 
motor control systems observed in healthy subjects and 
subjects with vestibular deficits [29]. When a more sophis- 
ticated representation is needed, the number of links can 
be increased to capture additional degrees of freedom. For 
example, Barin adopted multiple regression techniques 
and concluded that a two-link model is sufficient to ex- 
plain changes in center of pressure (COP) during postural 
control [32]. Kuo combined a two-link model with a full- 
state-feedback human sensorimotor control model and 
optimal control theory to study responses to small pertur- 
bations in the anterior-posterior (AP) direction [33]. Kooij 
et al. developed a three-link model based on optimal esti- 
mation theory to characterize the contribution of multi- 
sensory information to standing balance and concluded 
that a predictive element in the controller is essential to 
compensate for neural time delays [34]. Unfortunately, 
all of these models only capture movement in the AP dir- 
ection. Thus, Winter considered two separate models for 
movement in the AP and medial-lateral (ML) directions 
[35]; yet, even his models do not capture the dynamic coup- 
ling between the AP and ML motions, which is often sig- 
nificant during multidirectional perturbations. Finally, none 
of the models mentioned above incorporate biofeedback. 

To understand how biofeedback affects postural control, 
Goodworth et al. developed a vibrotactile biofeedback 
model for a single-link inverted pendulum model of stance 
[16]. In the single-link model, the ankle angle (which is 
controlled only through ankle torque) represents the com- 
bined kinematics of a number of body segments, i.e., the 
relative angle between the feet and legs as sensed by the 
proprioceptive system, the orientation of the head in space 
as sensed by the vestibular system, the relative position of 
the head with respect to the environment (which can be 
important in eyes-open scenarios although not considered 
in eyes-closed studies such as [16]), and the orientation of 
the torso in space as measured by the biofeedback device. 
Because single-degree-of-freedom models have only one 
sway angle to measure and only one joint torque to con- 
trol, representations of the sensory integration and con- 
trol processes associated with biofeedback are relatively 



straightforward. Although effective at capturing uni- 
directional postural responses to small perturbations, 
single-link models represent movement only along a 
single axis and are limited in their ability to capture 
multidirectional postural responses (especially during 
large perturbations). Furthermore, Goodworth et al.s 
representation of biofeedback is not readily scalable to 
models with higher degrees of freedom. 

This paper describes a new method to incorporate bio- 
feedback into a multi-degree-of-freedom model for human 
balance. Rather than considering an integration on the 
sensory side, the method integrates biofeedback after the 
existing postural control mechanism generates joint tor- 
ques and before they are actually applied to the joints. This 
scalable method is independent of the models degrees of 
freedom, how it handles sensory integration, and the 
modeling of its postural controller. The model is vali- 
dated against existing experimental data [28] to demon- 
strate its ability to replicate the experimentally observed 
average response trajectories of individuals with vestibu- 
lar deficits for three different perturbation directions 
and three different feedback device display configurations. 
The model is then used to predict the performance of a 
fourth display with a lower spatial resolution than those 
experimentally evaluated. Although this study features a 
vibrotactile feedback display, the unique approach pre- 
sented herein has the flexibility to describe auditory, 
visual, electrotactile, and multimodal feedback displays 
and scale to single or multi-segmented models of human 
stance, thereby rendering it suitable for a wide range of 
balance and balance rehabilitation applications. 

Methods 

The procedure followed in this study is outlined in Figure 1 
and explained in detail in this section. 

Mathematical modeling 
Biofeedback 

Typical biofeedback systems consist of a motion or force 
sensor to detect body kinematics or kinetics, respectively, 
a processor to estimate body kinematics or COP, and a 
feedback display to provide the user with an additional 
channel of information. This study assumed that the bio- 
feedback signal was either given in cardinal directions (i.e., 
in the AP and ML directions) by the device, or was 
decomposed by the individual into these cardinal direc- 
tions. Regarding the latter, if the biofeedback signal was 
given at an angle of 6 measured clockwise from the navel 
as viewed in the transverse plane, we assumed that the in- 
dividual decomposed the signal into its sagittal and coronal 
components according to cos 6 and sin 0, respectively. 

Accounting for experimental observations that indivi- 
duals control balance in the sagittal and coronal planes 
independently [18,39], the biofeedback signal was also 
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Figure 1 Procedural schematic of the methodology employed in this study. 
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assumed to be processed and utilized independently in the 
cardinal directions. Specifically, additional torque signals 
Tp and Tp were assumed to be generated due to bio- 
feedback independently for sagittal and coronal planes, 
respectively, as follows: 



K s 



T F S + 1 ' 



t£(s) 



T F S + 1 



(i) 



where Tp and Tp are the additional sagittal and coronal 
biofeedback torques, respectively, and are vectors with a 
dimension equal to the number of joints in the model. 
Thus, they are scalars in a single-link inverted pendulum 
representation of the body that models only the ankle joint, 
whereas they become three-dimensional vectors in a model 
that considers ankle, knee, and hip joints. The variable t f is 
the time constant associated with the reaction to the feed- 
back. It was assumed to be the same in both sagittal and 
coronal directions; hence the superscript S or C to differen- 
tiate the sagittal and coronal plane is omitted. The variable 
5 is the Laplace variable, a standard notation for represent- 
ing transfer functions in the frequency domain. Kp and Kp 
are the steady state magnitude vectors of sagittal and cor- 
onal torques due to feedback; i.e., Tf = K$ and Tp = Kp in 
the steady state. 

The biofeedback torques Tp and Tp were added to the 
joint torques generated by a postural controller as follows 



nrS 
1 PC 



Tf 



(2) 



where Tj represents the vector of joint torques and T PC 
represents the torque vector generated by the postural 
controller model. Figure 2 illustrates the integration of this 
biofeedback model into a generic feedback-based human 
postural control model. Note that the integration of bio- 
feedback in this approach is independent of the manner in 
which the human body, sensory integration, and postural 
control are modeled, and thus is easily applicable to multi- 
degree-of-freedom models without any constraints on the 
mathematical representations of the sensory and control 
systems. 

For validation purposes, the proposed modeling frame- 
work was applied to experimental results for a representa- 
tive biofeedback device as follows. The device used an 
inertial measurement unit to detect body motion and a 
vibrotactile display comprising 48 vibrating actuators (sub- 
sequently referred to as tactors) to provide a feedback sig- 
nal related to the measured torso tilt. The direction 6 of 
the feedback signal was measured clockwise from the 
navel and was equal to the torso tilt direction (azimuth), 
which was computed based on the arctangent of the AP 
and ML tilt components. The device feedback algorithm 
was such that the magnitude of the feedback signal was 
equal to the tilt angle (inclination) plus half the tilt angle 
rate-of-change (see [14] for details). Direction was displayed 
using as many as 16 columns of tactors based on a "nearest 
neighbor" principle in which the column closest to the ac- 
tual tilt direction is activated. The columns were placed 
around the torso at 22.5° intervals starting from the navel. 
Within each column, three rows of tactors encoded three 
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Figure 2 Block diagram representation of the biofeedback model. The variable s is the Laplace variable. 



different levels of feedback magnitude; no tactors were 
activated within a 1° dead zone and tactor activation 
progressed from the bottom to the top tactor row in a 
stepwise manner as sway increased. Tilt magnitude was 
modeled based on the algorithm used during the experi- 
mental studies (see [14] for details). K| and K§ were 
expressed as a function of the activated row r and the 
feedback direction 6 



( kf cosfl r = 1 
Kf(r,0) = < ldjcostf r = 2, 
I kf cosfl r = 3 



kf sinO 
1<2 sin# 
kf sin# 



(3) 



where k,, i = {1, 2, 3} for the three rows, are constant 
vectors. To ensure that biofeedback torque monotonic- 
ally increased with the activated row, the elements of k/ 
are constrained as /c 1>y < k 2 j < k 3 j. Here, refers to the 
; th element of the vector k t % 

Sienko et al. showed that this vibrotactile feedback device 
quickens subjects' return to upright and reduces sway fol- 
lowing a discrete surface perturbation [28]. Furthermore, 
they varied the number of active tactor columns to deter- 
mine the effect of spatial resolution on postural perform- 
ance [14,18,28]. For example, in a 3x4 (3-row, 4-column) 
display configuration, only four of the tactor columns 
(navel, spine, and left and right sides) were active, pro- 
viding 90° spatial resolution. No significant differences 



were observed among the display configurations, i.e., four 
columns were as effective as sixteen columns [28]. 

Bipedal stance 

The human body was modeled as a two-link inverted pen- 
dulum in three-dimensional space, where one link repre- 
sented the legs and the other link represented the upper 
body. The feet were not considered; rather, the body was 
assumed to be connected to the perturbation platform 
(described below) through the ankles. 

Standard vector second-order differential equations were 
used to express the rigid body dynamics of the two links 
[40]. Translational dynamics for each link were given by 
the following concise vector second-order equation 

F = mv (4) 

where F is the total external force acting on the link and v 
is the velocity of the center of mass of the link with respect 
to the inertial frame. Rotational dynamics were given by 



T = Id) + (o x la) 



(5) 



where T is the total external torque acting on the link, <o is 
the absolute angular velocity, and I is the moment of iner- 
tia matrix for the link coordinate frame at the center of 



mass. 



The ankles and hips were modeled as ideal spherical joints. 
No passive stiffness or damping effects were considered in 
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the joints, as the role of passive torque during stance has 
been reported to be negligible [41]. The knees were locked 
due to the fact that the perturbations were small This also 
helped reduce the degrees of freedom in the model 

The postural control mechanism was assumed to com- 
prise two independent linear state-feedback controllers for 
the AP and ML directions. The assumption of independ- 
ence is supported by experimental studies [18,39], and the 
linear state-feedback control assumption is standard in the 
literature [33,34]. State feedback was assumed to be with- 
out noise, but with delay. No feedforward or estimation 
mechanisms were considered. The equations for the pos- 
tural controllers had the same form in both the AP and 
ML directions, given by 



T PC (t) = -K x(t - t d ) 



(6) 



where T PC (t) = [T PCa (t) T PC h(t)} T is the vector of postural 
control torques at the ankle and hip, K is a 2x4 matrix of gains, 
andx(£- t d ) = \a a (t- t d ) a h (t - t d ) a a (t - t d ) a h (t - t d )] T 
is the state vector of ankle and hip angles and angular vel- 
ocities with a time delay of 

Sensory integration was neglected for simplicity. The 
states from the human body model were directly fed into 
the postural controller; i.e., the Sensory Integration block 
in Figure 2 was not used. 

The model was implemented in 20-Sim (Controllab 
Products B.V., Enschede, the Netherlands). Inputs to the 
model were the discrete surface perturbations used in 
[28], which are described next. 

Discrete surface perturbations 

Perturbations were applied using a custom-built 2.1 m 2 
platform [42,43] that could move in an earth-horizontal 



plane by independent control of motion in two orthog- 
onal (x and y) directions. Figure 3 shows the platform 
velocity inputs for the three perturbation directions con- 
sidered in this study (90°, 180°, and 225°). Perturbations 
comprised a constant acceleration phase for 100 ms, fol- 
lowed by 200 ms constant velocity and 100 ms constant 
deceleration phases. The perturbation magnitude was sub- 
ject specific and ranged from 50 to 70 mm. Two-axis tilt 
(roll and pitch), COP, and platform position were collected 
at 100 Hz. Additional details about the perturbations can 
be found in [28]. The average control inputs to the plat- 
form from the experimental study were used in the simu- 
lations described below. 

Model parameterization 

The model was parameterized using experimental and an- 
thropometric data from six subjects with bilateral vestibu- 
lar hypofunction who participated in a study aimed at 
characterizing the effect of multi-directional vibrotactile 
biofeedback on postural stability during discrete multidir- 
ectional support surface perturbations [28] . Detailed infor- 
mation about the subjects can be found in [14,28]. 

The average subject height and mass was 1.78 m±0.09 m 
(SD) and 86 kg ±5.4 kg (SD), respectively. Body link mo- 
ments of inertia were calculated based on standard formu- 
las for the cylinders and ellipsoids used to approximate the 
shape of the subjects' legs, torsos, and heads. The masses 
and lengths of the body links were found using percent-of- 
total-weight [44] and percent-of-total-height [45] data. 

Postural controller parameters were found by minimizing 
the integral of the square of the sum of the normalized dif- 
ferences between the experimental and simulated trajector- 
ies of the AP and ML sway and COP data. Normalization 



E 

5 



Time (s) 



225 c 



180° 




co 0-2 
E 



Time (s) 

-*-90° 



4 

Time (s) 



Figure 3 Perturbations. Bird's-eye view showing the perturbation directions, the positive directions for the x- and y-axes, and the platform 
velocity inputs for each perturbation direction. 
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was performed using the corresponding peak experimental 
quantities. Hence, the optimizations sought to minimize 
the following function 
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where 0 represents sway, COP is the center of pressure, 
subscripts exp and sim denote experimental and simula- 
tion data, respectively, and superscripts S and C denote 
the sagittal and coronal planes. AP sway was not included 
in Eq. (7) for 90° perturbations, and ML sway was not 
included for 180° perturbations, because the model exhi- 
bits a unidirectional response in these cases. Different 
control gains were used for the three perturbation direc- 
tions to ensure best fit. 

The bipedal stance model parameters were held con- 
stant while tuning the biofeedback model parameters for 
the 4-, 8-, and 16-column display models using the averaged 



experimental data from the 4-, 8-, and 16-column dis- 
play configuration studies, respectively, and minimizing the 
same objective function given in Eq. (7). Only the 225° per- 
turbation direction was considered for the 8- and 16- 
column displays, because they are identical to the 4-column 
display for the 90° and 180° perturbation directions. 

A 1x2 configuration was also considered to determine 
if a configuration simpler than those tested experimen- 
tally could be as effective, recognizing that such a con- 
figuration could correlate with fewer electromechanical 
components and lower device cost. In this configuration, 
the two columns were assumed to be aligned with the AP 
axis. The parameters for this configuration were tuned 
using the 3x4 experimental data. The 90° perturbation dir- 
ection was not considered, since the 1x2 configuration 
does not provide any feedback in the ML direction. 

Minimizations were done by a parameter sweep in the 
design space followed by using the best point as an initial 
guess for an optimization with the Perpendicular Search 
approach. In cases where the design space was large, such 
as the 225° perturbations due to the simultaneous tuning 
in both sagittal and coronal controllers, a sensitivity ana- 
lysis was performed to reduce the parameter sweep space 
for tractable computation. However, all of the parameters 
were included in the subsequent optimization. 
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Model validation metrics 

The following metrics were used to quantify the models 
quality of fit: 

1. Average error: absolute value of the instantaneous 
difference between the experimental response and 
the simulated response averaged over 3 s after 
perturbation; values near zero indicate a good fit. 

2. Cross-correlation value: cross-correlation value 
between the experimental response and the 
simulated response over 3 s after perturbation; values 
near 1.0 indicate a good fit. 

The duration of the time window was chosen as 3 s to 
capture the entire transient response following perturb- 
ation. The transient part of the response is of interest, 
because in the steady state the model returns to the ini- 
tial states of zero sway and COP. 

Results 

Figure 4 shows the experimental and simulated trajec- 
tories for the tactors-off case. Figures 5 and 6 show the tra- 
jectories for the 3x4, 3x8, and 3x16 display configurations. 
Figure 7 shows the simulated 1x2 display configuration and 



3x4 experimental trajectories for the 180° and 225° perturb- 
ation directions. The average error and cross-correlation 
values are summarized in Table 1. 

While obtaining these fits, the maximum ankle and hip 
torques were observed to be approximately 75 N-m and 
25 N-m, respectively. As an example, joint torque trajec- 
tories for the 225° perturbation direction are shown in 
Figure 8 for the tactors-off and 3x4 configurations. Dur- 
ing the tactors-on conditions, the maximum additional 
ankle and hip torques due to biofeedback were about 
8 N-m and 3 N-m, respectively. It was observed that no 
additional coronal torques were needed to obtain the 
best fits in the tactors-on conditions, including the 90° 
perturbations. The joint torques obtained by fits to the 
model were biomechanically feasible [46-48] and con- 
sistent with previously reported values [49] . 

Based on these results, the bipedal stance model com- 
prising a multibody model and two independent full- 
state-feedback controllers for the AP and ML directions 
was considered valid for representing human stance 
during discrete multidirectional support surface pertur- 
bations, and the biofeedback model was considered 
valid for capturing the effect of vibrotactile biofeedback 
on balance. 
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Figure 5 Average experimental and simulated responses for the 3x4 display configuration. 
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Discussion 

Sensory reweighting [11] and sensory addition [16] have 
previously been described as potential mechanisms by 
which biofeedback devices supplement native sensory 
inputs to inform corrective motor torques and thereby 
decrease postural sway. Sensory reweighting is a general 
term used to describe both the real-time (i.e., while bio- 
feedback device is being worn during a single training 
session) and long-term/plastic (i.e., post extensive train- 
ing in the absence of wearing the biofeedback device) 
effects of supplemental information on the gains of native 



sensory inputs. In other words, the improved quality of 
the information about the body's motion with respect to 
the gravito-inertial vector supplied by the non-native 
supplemental channel of information would result in an 
increased reliance on native inputs that correlate highly 
with this information. This approach infers upstream in- 
tegration of the supplemental information within the 
central nervous system (CNS) (e.g., non-native informa- 
tion from the biofeedback device may be integrated with 
native inputs upstream within the multimodal sensory 
association areas in the cerebral cortex). However, to date, 
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Table 1 Model validation metrics (average error and 
cross-correlation values) across display configurations 
and perturbations 

Average error (Sway) 



AP ML 





180° 


225° 


225° 


90° 


Tactors off 


0.16° 


0.19° 


0.09° 


0.19° 


1x2 


0.85° 


0.27° 


0.11° 




3x4 


0.33° 


0.36° 


0.10° 


0.18° 


3x8 




0.21° 


0.14° 




3x16 




0.30° 


0.12° 








Cross-correlation value (Sway) 








AP 


ML 






180° 


225° 


225° 


90° 


Tactors off 


0.998 


0.99 


0.97 


0.95 


1x2 


0.94 


0.97 


0.98 




3x4 


0.98 


0.95 


0.98 


0.96 


3x8 




0.98 


0.95 




3x16 




0.97 


0.93 








Average error (COP) 








AP 


ML 






180° 


225° 


225° 


90° 


Tactors off 


0.47 cm 


0.21 cm 


0.33 cm 


0.46 cm 


1x2 


0.49 cm 


0.26 cm 


0.51 cm 




3x4 


0.46 cm 


0.23 cm 


0.51 cm 


0.42 cm 


3x8 




0.34 cm 


0.49 cm 




3x16 




0.27 cm 


0.42 cm 








Cross-correlation value (COP) 








AP 


ML 






180° 


225° 


225° 


90° 


Tactors off 


0.98 


0.99 


0.98 


0.98 


1x2 


0.97 


0.98 


0.93 




3x4 


0.96 


0.99 


0.93 


0.98 


3x8 




0.97 


0.92 




3x16 




0.98 


0.95 





Results for all display configurations, with the exception of the 1x2 case, were 
obtained by fitting the model to the corresponding experimental data. The 
model for the 1x2 display configuration was fit to the experimental data from 
the 3x4 case. All parameters obtained for the tactors-off condition were held 
fixed in the tactors-on configurations. 



there is limited or no evidence in the literature to support 
either the real-time or long-term/plastic sensory reweight- 
ing scenarios. 

Sensory addition is based on the notion that the infor- 
mation from the non-native channel is added to the infor- 
mation from native channels, with the gains of the native 
channels remaining unchanged. The modeling results 
reported by Goodworth et al. [16] suggest that the infor- 
mation obtained from their vibrotactile biofeedback device 



added to the native vestibular and proprioceptive feedback 
without changing the reliance on the native sensory inputs 
(i.e., the vestibular and proprioceptive gains remained un- 
changed). However, it should be noted that this approach 
also requires upstream integration of the supplemental in- 
formation with the native inputs within the CNS. 

This paper contributes to the literature by integrating 
biofeedback into a postural control model as an additive 
motor torque. Our approach integrates biofeedback on 
the torque side, unlike sensory reweighting or addition 
approaches that integrate biofeedback on the sensory side 
(i.e., within the Sensory Integration block in Figure 2). 
Our approach is important for three reasons. First, there 
is a lack of evidence for either real-time or long-term sen- 
sory reweighting. Second, the device used in this and our 
previous studies [7,11,14,16,18], which requires subjects to 
produce volitional (cognitively processed) responses to the 
supplemental information provided, suggests that it is more 
realistic to model the addition of the non-native body mo- 
tion cues further downstream within the CNS, such as 
within the motor association cortex or primary motor cor- 
tex. Third, this approach is compatible with single or mul- 
tiple degree-of-freedom representations of the body, 
multiple sensory integration models [16,29,33,34] (although 
not explored in this study), and various feedback modal- 
ities. Although both the sensory addition biofeedback 
model of Goodworth et al. and the additive torque bio- 
feedback model described herein fit their respective ex- 
perimental data well, we believe that the flexibility and 
scalability of the additive torque biofeedback model make 
it more suitable for a wide range of applications. 

While looking for the simplest model to represent 
human stance in this study, we also considered a single- 
link representation. However, even though this repre- 
sentation was reported to be successful by Goodworth 
et al. [16], we found that it failed to adequately capture the 
experimental responses we considered. Specifically, the AP 
sway average errors were about nine times larger than those 
obtained with the two-link model. We postulate that the dif- 
ference in adequacy of the single-link representation may be 
due to the difference in the types of perturbations. In par- 
ticular, Goodworth et al. used continuous rotary platform 
perturbations, whereas we obtained our experimental data 
using discrete translational support surface perturbations. 
Hence, a direct comparison with the sensory addition 
model was not readily possible. 

The two-link biomechanical model with associated 
controllers and biofeedback scheme is capable of captur- 
ing experimental response trajectories with low average 
errors and high cross-correlation values in both the AP 
and ML directions for all perturbation directions and de- 
vice display spatial resolutions. Since the optimal fits to 
experimental data were obtained without the need for 
any coronal biofeedback torques, we conclude that the 
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Figure 8 Simulated joint torque trajectories for the 225° perturbation. 



3x4 factor configuration 




benefit of increasing the display resolution is question- 
able. Thus, our simulations support the experimental 
findings [14,28] that higher resolution displays do not ne- 
cessarily correspond to better performance, possibly due 
to: 1) the larger effective stiffness in the ML direction than 
in the AP direction; and 2) the feedback display scheme. 
Specifically, the multidirectional perturbations elicited 
greater body movement in the AP than in the ML direc- 
tion. This phenomenon, coupled with the "nearest neigh- 
bor" display scheme used in this paper, means that 
displayed information typically aligns with the AP direc- 
tion for all display resolutions. Thus, the increased reso- 
lution has a negligible effect. 

In fact, the simulation study with the 1x2 configuration 
showed that this configuration can match the experimen- 
tal data for the 225° perturbation well and hence can be as 
effective as the 3x4 configuration for the 225° perturb- 
ation. However, for the 180° perturbation, the body sway 
fit obtained for the 1x2 configuration was not as good as 
the other fits (see Table 1), possibly because the 180° 
perturbation causes a larger body sway amplitudes. 
Therefore, a higher number of rows could be beneficial 
for larger sway amplitudes. These results should encour- 
age further studies with lower display resolutions [50]. 
In the simulations for the 3x4, 3x8, and 3x16 configura- 
tions we observed that no biofeedback torque was gen- 
erated in the hips when the first row of the device was 
active; instead, it caused biofeedback torques only in the 
ankles, whereas activation of the second row caused 
biofeedback torques in the hips. These observations align 
with previous observations that standing human subjects 
employ the "ankle strategy" for small perturbations and 
the "hip strategy" for larger perturbations [33]. 

In the simulations for the 3x4, 3x8, and 3x16 config- 
urations, we observed that activation of the first row of 
the device caused biofeedback torques only in the ankles, 
whereas activation in the second row caused biofeedback 
torques also in the hips. These observations align with 
previous observations that standing human subjects 
employ the "ankle strategy" for small perturbations and 
the "hip strategy" for larger perturbations [33]. 



The fact that the model can account for the experi- 
mental tactors-on responses without changing the para- 
meters of the full-state-feedback controllers of posture 
concurs with previous findings that biofeedback does not 
necessarily have to act through sensory substitution, or 
cause sensory reweighting to have an impact on real-time 
balance performance [16]. 

This paper s limitations are as follows. A strict focus on 
perturbed stance indicates that the model and results may 
not be immediately transferable to studies of gait. Our 
model is deterministic in nature and thus does not capture 
the variability of the response in its current form. Because 
the model was developed and validated for small perturba- 
tions that do not elicit a need to move the arms or the feet 
and do not elicit nonlinear postural responses, studies 
with more severe perturbations will likely require more 
complicated models with nonlinear control schemes. Fi- 
nally, the model does not consider the cognitive load asso- 
ciated with using the feedback device. 

We believe that our biofeedback modeling approach has 
broad applicability within the field of balance rehabilitation. 
While this study considered a specific (vibrotactile) bio- 
feedback device, the model is not limited to this feedback 
modality. Our model is easily modified to capture postural 
responses to auditory, visual, electrotactile, or multimodal 
feedback displays by incorporating different reaction time 
constants (t f ) and feedback gains (K F ). Although we have 
illustrated the models utility for studying display reso- 
lution, we suggest that the approach can be used to study 
different device feedback algorithms and factor activation 
schemes. The application of the model can also be 
extended to a wider range of subject groups, e.g., preg- 
nant, obese, and aging populations. Such studies are be- 
yond the scope of this paper, but are important directions 
for future research. 

Conclusions 

The effect of biofeedback on body sway trajectories during 
perturbed stance can be accounted for by modeling the 
biofeedback as a torque signal that is added to the existing 
joint torques generated by the postural controller. Unlike 
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the sensory addition approach, this torque addition ap- 
proach is independent of the way the body, sensory integra- 
tion, and postural control are modeled and hence provides 
a scalable method to integrate biofeedback. The proposed 
model suggests that biofeedback can work without neces- 
sarily requiring a sensory reweighting or substitution. For 
the specific validation study performed, the model also sug- 
gests that increased resolution in vibrotactile biofeedback 
displays does not necessarily lead to better performance in 
terms of reduced body sway. In fact, providing feedback in 
the sagittal plane may be adequate even during small multi- 
directional perturbations. 
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